Nodal Liquid Theory of the Pseudo-Gap Phase of High- Tc Superconductors 
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We introduce and study the nodal liquid, a novel zero-temperature quantum phase obtained by 
quantum-disordering a d-wave superconductor. It has numerous remarkable properties which lead 
us to suggest it as an explanation of the pseudo-gap state in underdoped high-temperature super- 
conductors. In the absence of impurities, these include power-law magnetic order, a T-linear spin 
susceptibility, non-trivial thermal conductivity, and two- and one-particle charge gaps, the latter evi- 
denced, e.g. in transport and electron photoemission (which exhibits pronounced fourfold anisotropy 
inherited from the d-wave quasiparticles) . We use a. 2 + 1-dimensional duality transformation to de- 
rive an effective field theory for this phase. The theory is comprised of gapless neutral Dirac particles 
living at the former d-wave nodes, weakly coupled to the fluctuating gauge fleld of a dual Ginzburg- 
Landau theory. The nodal liquid interpolates naturally between the d-wave superconductor and the 
insulating antiferromagnet, and our effective field theory is powerful enough to permit a detailed 
analysis of a panoply of interesting phenomena, including charge ordering, antiferromagnetism, and 
d-wave superconductivity. We also discuss the zero-temperature quantum phase transitions which 
separate the nodal liquid from various ordered phases. 



I. INTRODUCTION 



The discovery oLthe cuprate high-temperature super- 
conductors in 1986EI was a watershed in the recent history 
of condensed matter physics, an event which stimulated 
intense experimental and theoretical activity. As sample 
quality and experimental precision have advanced, these 
materials' rich jphase diagram and phenomenology have 
come into focus.0 However, many theoretical efforts have 
not reached fruition because of a serious obstacle, namely, 
that these materials are apparently in a strongly-coupled, 
non-perturbative regime. To put it more bluntly: there 
is no obvious small parameter which facilitates an expan- 
sion about a tractable model, fn this paper, we promul- 
gate the existence of a weakly-coupled 'dual' description 
of a particularly exotic region of the phase diagram. This 
description paves the way for controlled calculations of 
experimentally measurable quantities. 

Before plunging into our exegesis of this dual descrip- 
tion and its consequences, let us briefly review the phe- 
nomena which we wish to explain. The phase diagram as 
a function of temperature T and doping x (the nature of 
the doping varies from material to material, but is gener- 
ally believed to be proportional to the effective hole con- 
centration in each Cu02 layer) is indicated schematically 
in Fig. 1. Best understood and in many cases very well 
characterized are the undoped {x = 0) materials, which 
are insulating antiferromagnets (AFs) below the Neel 
temperature T/v. At moderate dopings (0.1 ^ x ^ 0.3) 
and low temperatures (T < Tc{x)), superconductivity 
occurs. For many phenomenological purposes the super- 
conducting phase is adequately described by the same 
Ginzburg-Landau and London theory standard for con- 
ventional superconductors (SCs).B One important dis- 



tinction between the high-Tc and conventional SCs, sus- 
pected for many years and now generally accepted i^er 
a number of beautiful and compelling experiments, uu is 
their d^^-y^ (d-wave) pairing symmetry. This pairing 
symmetry is a crucial ingredient in a zero-temperature 
quantum description. In particular, d-wave symmetry 
leads to gapless quasiparticles residing at the four nodes 
of the pair wavefunction. 



T*(x) 



Pseudo-Gap 




Fig. 1: Schematic phase diagram of a high-temperature 
superconductor as a function of doping x and tempera- 
ture T. 



Recent studies have revealed puzzling behavior 
in the underdoped region between the AF and 
SC. Below the high-temperature dotted line {T*{ 
in Fig. Ip. anglejjpiolved photoemissiorj_.(pAJlPES) 
transport ,eI NMR,El'E2l optical conductivityjlJ'E3 and other 
measurementffl indicate a dramatic reduction of low- 
energy (single-)electronic and spin degrees of freedom. 
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Furthermore, the ability of ARPES measurements to re- 
solve wavevector dependence exposes an angular varia- 
tion similar to that of d-wave quasiparticles in the SCing 
phasaj'tD. This portion of the phase diagram is commonly 
called the pseudo-gap region. The ultimate nature of the 
corresponding underlying quantum ground state is an in- 
triguing theoretical puzzle, and a principal subject of this 
paper. 

To proceed, we look to the experiments for guidance. 
They indicate three types of excitations which are im- 
portant below the dotted line in Fig. 1: the ordering 
fields related to antiferromagnetism and superconductiv- 
ity, and d-wave quasiparticles near the four nodes. Con- 
spicuously absent from this list are electrons and holes at 
an ordinary Fermi-surface. The physics of this omission 
is that pairing occurs (due to unspecified strong interac- 
tion physics) at the high energy T*{x). Given these in- 
gredients, one natural theoretical strategy is to attempt 
to approach the pseudo-gap state by increasing x from 
the AF at half-filling. Many researchers have already 
attempted this approach, but it remains inconclusive. 

We, instead, tackle the pseudo-gap state from the 
right, literally. To do this, we must contemplate 
quantum-disordering the d-wave superconductor. For 
simplicity, we will assume for the moment a purely two- 
dimensional model of a single Cu02 plane. We imagine 
that pairing establishes a local superconducting d-wave 
order parameter A(x, t) = |A|e"'', where x is the two- 
dimensional coordinate and t is time. The experimen- 
tal properties of the pseudo-gap state can be interpreted 
as an indication that |A| is large in this region, so that 
quantum fluctuations of the phase of the order-parameter 
tf must be responsible for the lack of off-diagonal long- 
range-order ((A) = 0), even as T ^ 0. The important 
long-distance dephasing is accomplished by vortex loops 
and lines, around which Lp winds by ±27r. To destroy 
the long-range correlations in A, we must unbind vortex 
loops of arbitrarily large size, just as vortcx-antivortex 
pairs unbind above the Kosterlitz-Thouless transition 
temperature in a two-dimensional superfiuid. To imple- 
ment this unbinding, it is axt remely helpful to use the 
2 -I- 1-dimensional duaiitjt^lllfl relating an XY-model for 
A and a Ginzburg-Landau model with complex field 
("disorder parameter") interacting with a gauge field a. 
The duality interchanges Cooper pairs and vortices, so 
that the desired quantum disordered d-wave state is the 
ordered (condensed) phase, ($) ^ 0, of the Ginzburg- 
Landau theory. The imaginary-time effective action for 
this dual theory is nothing but the Ginzburg-Landau free 
energy functional for a three-dimensional superconductor 
at finite-temperature. 

The fate of the d-wave quasiparticles in this construc- 
tion requires particular care, since these ar e stro ngly cou- 
pled to the fluctuating pair field. In Sec. II A, we show 



spin-ful) quanta we call nodons. In the dual variables, 
the fundamental Lagrangian of our model is 

C = ^j\[idt — vpT^idx — v^T^idylipi + (1 <-» 2, a; ^ y) 

+2XJo + Cn. (1.1) 

In the sections which follow, we will elucidate the physics 
of this Lagrangian in some detail, so we restrict our- 
selves, in this introduction, to a whirlwind tour. The 
fields ipjjipj, j = 1,2 are the nodon creation and anni- 
hilation operators at the two antipodal pairs of nodes. 
J\ is a bilinear in the nodon operators which has an in- 
terpretation in the d-wave superconducting phase as the 
quasiparticle electrical 3-current. (f> is the complex scalar 
field representing the vortices, and is the gauge field 
which is dual to the phase of the superconducting or- 
der parameter. The term proportional to A describes the 
effects of particle/hole asymmetries, and the Cn term de- 
scribes the coupling of the nodons to antiferromagnetism, 
which we will return to presently. 

A remarkable result of calculations with Eq. 



1.1 



is that 

gapless nodons survive the quantum disordering of the 
SC! The nodons are like the smile of the Cheshire cat: 
the dx2_y2 order parameter is gone, but the nodes re- 
mai n. T he consequent Nodal Liquid (NL) described by 
Eq. 
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how the quasiparticles can be treated by extracting the 
C7(l) phase from the bare electron operators. Once the 
phase is extracted, one is left with a set of gauge-invariant 
fermion operators which create electrically neutral (but 



is a distinct and novel zero-temperature quan- 
tum phase with a number of fascinating properties. For 
simplicity, consider first a hypothetical NL phase at half- 
filling in the absence of impurities. The possibility of AF 
ordering will be included later via Cn- Wc hypothesize 
that antiferromagnetism might be avoided and the NL 
achieved in a half-filled square lattice model by adjusting 
an attractive nearest-neighbor interaction and second- 
neighbor electron hopping amplitude. The NL is a nomi- 
nally insulating state, with non-zero gaps Ai and A2 for 
adding both individual electrons/holes and Cooper pairs, 
respectively. Gapless nodons with anisotropic ballistic 
dispersion (lu ^ k), however, persist, and can carry both 
spin a nd c urrent. With particle/hole symmetry (A = 
in Eq. IT), we expect power law (~ l/la;!"*) spin cor- 
relations at (0,0), (tTjTt), (77,0), and symmetry-related 
points in the Brillouin zone. Scaling arguments lead us 
to expect a weak dissipative dynamic contribution to the 
conductivity which, in the presence of particle/hole sym- 
metry, would vary as Re(j{uj,T = 0) ~ uj^. Nodons 
also contribute a quadratic specific heat Cnodon ~ aT^. 
Despite the similarity of the nodons to d-wave quasipar- 
ticles, the single-particle spectral function is predicted to 
show a gap at the Fermi energy in the NL. We expect, 
however, this gap to be strongly angle-dependent: of or- 
der the pairing scale T* away from (±7r/2, ±7r/2) and 
reduced to Ai ^ T* at these special points. A parti- 
cle/hole asymmetric NL should exhibit similar behavior, 
but with singularities shifted from (±7r/2, ±7r/2) in mo- 
mentum space, and an even smaller contribution to the 
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low-frequency dissipative conductivity which, in the ab- 
sence of umklapp scattering and impurities must come 
from interactions with phonons. 

Consider next doping the NL. Naively, this can be 
modeled via an increase in the in-plane hole concentra- 
tion, although the actual transfer of charge to the Cu02 
layers may be not be complete. In the grand canonical 
ensemble, charge is added by increasing fj, = Aq/2. In 
the dual theory (see Eq. 1.1), this chemical potential acts 
like an external magnetic field /i = in the Ginzburg- 
Landau theory (but, unlike a magnetic field, it of course 
does not break time-reversal invariance). For small /i, the 
system remains in the Meissner phase, and no dual flux 
penetrates - i.e. no charge is added to the system within 
the charge gapm Following the analogy with Ginzburg- 
Landau theory,cl we expect that the nature of doping de- 
pends upon the Abrikosov parameter = ^v/^v^ where 
and are effective dual penetration and coherence 
lengths, respectively. For Ky > l/"\/2, type II doping 
occurs, and the "field" penetrates first for /i > ^ci in 
a dual flux lattice. Dual flux tubes are in fact Cooper 
pairs, so this is a paired Wigner crystal (PWC) state. 
We show in Sec. Ill that gapless nodons survive the dop- 
ing and coexist with the PWC. We generally expect the 
displacements of the Wigner crystal to be pinned either 
by the periodic lattice potential or disorder (either is ef- 
fective when arbitrarily weak), so that this phase remains 
insulating. As doping increases from zero, the character- 
istic nodon momenta shift further from (±7r/2, ±7r/2). 
Continued doping to /Lt > /Zc2 leads to another transi- 
tion into the "normal" state of the dual theory, which is 
nothing but the d-wave SC. Neglecting disorder, and with 
weak lattice effects, quantum fluctuations are expected to 
drive this 2 -I- 1-dimensional flux lattice melting transition 
weakly first order. For k ^ l/v^, one has instead type I 
doping at a single "critical field" /ic- This is a first order 
transition, accompanied by a jump in the hole concentra- 
tion from to Xc at fi = ^c- In the canonical ensemble 
with fixed < x < Xc, one expects two-phase coexis- 
tence, i.e a "mixed" phase in the dual Ginzburg-Landau 
theory. Taking into account long-range Coulomb inter- 
actions, one arrives at the frustrated phase sepaxation 
physics discussed at length by Emery and Kivelson,Ej and 
all the consequent issues. Crude arguments (see Sec. HI) 
suggest K ~ 1 — 3, in the weakly type II limit, but close 
enough to the threshhold value to allow for different sce- 
narios in different materials. Regardless, type II and 
some type I schemes imply NL phases at finite doping 
before the onset of superconductivity. 

An intermediate NL phase is extremely appealing from 
the theoretical point of view, as it offers a compelling 
interpolation between the undoped AF and the d-wave 
superconductor. Consider the following three important 
energies : the single-particle gap, Ai, the minimum en- 
ergy required to add a charge ±e and spin s = 1/2 to the 
system; the two-particle gap, A2, the minimum energy 
to add charge ±2e and spin s — 0: and the spin gap, A^, 
the energy required to add spin s = 1 but no charge. In 



the AF, both single- and two-particle gaps are non-zero, 
but the spin gap vanishes due to low frequency magnons. 
In the d-wave superconductor, the two-particle gap van- 
ishes, since the pairs have condensed, but spin and single- 
particle gaps are "almost" non-zero (which we call "O'^"), 
since only the quasiparticles carry these quantum num- 
bers, and their density of states vanishes with the energy. 
Passing from the SC to the NL, A2 changes from zero to 
non-zero, and Ai changes from 0^ to a true gap. The 
transition to the AF occurs then simply by developing a 
non-zero staggered magnetization. 

This transition and other magnetic physics is discussed 
in Sec. IV, using the Neel Lagrangian density 



CN = ^\d,N\^-VM{\N\) 

+,g[N- V'VV^'cTTAt +h.c.]. 



(1.2) 



where g measures the strength of the coupling between 
the Neel order parameter and the nodons. Eq. |1.2| can 
be obtained by introducing a 2fci? density-density interac- 
tions between the nodons and decoupling the antipodal 
terms with the Neel vector N. Let us once again con- 
sider first the case of half- filling with particle/hole sym- 
metry. For sufficiently strong interaction g, or when the 
quadratic coefficient in Vn^N) is negative, one obtains 
an AF phase with (N) 7^ 0. In this phase the nodons 
develop a gap and low-energy spin quanta are carried en- 
tirely by spin waves. Depending upon the "mass" r$ of 
V$ , this is either a simple AF or AF order coexisting with 
a d-wave SC. Decreasing g or increasing vn destroys the 
long-range AF order and liberates the nodons. This inter- 
esting phase transition is discussed in Sec. IV. Increasing 



r$ results in a further transition to the d-wave super- 
conductor, which we believe is in the three-dimensional 
inverted- XY universality class. Tuning r$ = r^r = de- 
scribes a multicritical point connecting directly the AF 
and d-wave SC phases. 

Without particle/hole symmetry, A 7^ 0, and another 
possible phase exists: the coexisting AF and Nodal 
Liquid (AF/NL), with long-range AF order at (7r,7r) 
and power-law spin-density-wave correlations from the 
nodons at incommensurate wavevectors. This phase may 
be difhcult to distinguish experimentally from the pure 
AF, and it seems possible that some of the well-studied 
undoped cuprate materials might well be in the AF/NL 
phas e rather than the pure AF. In any event, the model of 
Eqs. 



1.1 



1.2 



provides a simple basis for understanding the 
suppression of Neel order upon doping. To see this, as- 
sume that at half-filling rjy > 0. Agreement with experi- 
ment then requires that g be sufficiently strong to induce 
AF order. As the chemical potential fj, is increased above 
the charge gap to induce holes into the system, the hole 
density or dual "internal field" becomes non-zero. From 
Eq. 1.1, this creates an effectively larger particle/hole 
asymmetry Aoff — A ~ a;. This presents a competition. 
By ordering the Neel vector the system can create a gap 
for the nodons and reduce their kinetic energy. However, 
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at finite Aeff the nodon Dirac point would prefer to move 
away from (±7r/2,±7r/2), which reduces the gain in ki- 
netic energy. As x increases, therefore, we may expect 
to drive transitions from the AF to AF/NL and pure 
NL phases. Of course, there are in fact many different 
scenarios for type I and type II doping, smaU or large in- 
trinsic A, etc. These are discussed in Sec. IV. Once mag- 
netism has been discussed, we conclude with a summary 
of the main points of the paper, open issues, relations to 
other work, and a brief discussion of experimental impli- 
cations. Finally, two appendices include technical details 
of microscopic and e-expansion calculations. 



II. D-WAVE AND DUALITY 

A. Model and Symmetries 

Consider a tight binding model of electrons hopping 
on a square lattice, with a local Hamiltonian satisfying 
certain general symmetries. We will assume the system is 
both U{1) and SU{2) invariant (i.e. we neglect spin-orbit 
coupling), and has time- reversal, reflection, and four- fold 
rotational symmetry. Sometimes it will also be conve- 
nient to specialize to models which possess an additional 
discrete particle/hole symmetry. We denote lattice elec- 
tron creation and annihilation operators as c^(a?) and 
Cq,(x), where x is the two-dimensional coordinate in the 
frame with x — xi and y — X2 parallel to the a and b 
crystalline axes (i.e. to the Cu~0 bonds). Here a is a 
spin label. In momentum space the kinetic energy takes 
the usual form, 



Ho 



kot 



(2.1) 



and at this stage we allow for general electron interac- 
tions: 



Hi, 



E 

k,q,k' 



V{k, q, k')cl^^^^^clpC-k'+q,0Ck',a- (2.2) 



by 



A discrete particle/hole transformation is implemented 



(2.3) 



with 7? = (tt, tt). Many common models (e.g. the Hub- 
bard and t-J) are invariant under a particle/hole trans- 
formation at half-filling. Invariance of the kinetic energy 
implies that efc — —ek+s, a form valid with near neigh- 
bor hopping. However, a second neighbor hopping term 
violates particle/hole symmetry. 

As discussed in the introduction, we wish to describe 
the physics below the relatively strong d-wave pairing 
scale T* , in order to approach the pseudo-gap phase from 
the superconducting side. To do so, we imagine introduc- 
ing a d-wave order parameter 



^kiq) = X! ^(^' 1' k'){C-k'+qlCk'^)- 



(2.4) 



BCS theoryES can be implemented in terms of the spa- 
tially varying pair field, obtained via Fourier transfor- 
mation, Ak{x). The self-consistent gap equation is usu- 
ally solved for a spatially uniform order parameter, with 
Afc = Ak{x). Singlet pairing implies Afc — A_fc, and in a 
d-wave superconductor has four zero's or nodes as k 
varies around the Fermi surface. Our strategy will be to 
obtain an effective field theory which has a local d— wave 
gap function, determined by strong coupling physics be- 
low some length scale A~^ of say 5-10 lattice spacings, 
but which can fluctuate quantum mechanically on longer 
spatial scales. These longer length scale quantum fluc- 
tuations will be responsible for quantum disordering the 
d-wave superconductor, and will allow us to access a new 
phase - the nodal liquid. As we shall see, an important 
role is played by the d— wave quasiparticles, which survive 
the quantum disordering. To implement this approach, 
we first briefly recapitulate the properties of quasiparti- 
cles in the d— wave superconductor. 



B. Quasiparticles 

With spatially uniform d-wave order given by A^, 
the effective Hamiltonian for the quasiparticles is H — 
Hq + Hi , with Hq the kinetic energy and 



Hi=J2[^kclcl,^+A 



(2.5) 



Since A^ = A^^k for singlet pairing, it is natural to break 
sums into positive and negative ky. To do so, consider a 
four component fermion field, Taa{k), at each wavevector 
with ky > positive: 



"^aaik) — 



Til 
T21 
T12 
'T22 



-fei 



(2.6) 



In the second column vector the minus sign has been 
introduced so that Taa transforms like a spinor un- 
der an SU(2) rotation, i.e. as Taa — > UapTafS- Here 
U = exp(z0 • cr) is a global spin rotation with Pauli ma- 
trices (Tap- 

In these variables, the quasiparticle Hamiltonian be- 
comes 

Hqp - 5]' Tt(^)[r-efe + r+A^ + r-A^]T(fc), (2.7) 



where the prime on the summation denotes over ky pos- 
itive, only, and we have introduced a vector of Pauli ma- 
trices. Tab acting in the particle/hole subspace. Also, we 
are employing the notation — (r^ ± ir*')/2. 
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and the particle/hole symmetry breaking term, 



■^2 



Fig. 2: The wavevectors Ki, q^, Qy in relation to the 
a,b axes. The dotted hne represents the putative Fermi 
surface. 

With approximate particle/hole symmetry, the d-wave 
nodes are located near the special wavevectors iKj , with 
Ki — (7r/2, 7r/2) and K2 = (— 7r/2, 7r/2). Since our aim 
is to obtain an effective description at low energies and 
long lengthscales, it is sufficient to focus on the gapless 
modes near these points, integrating out the electrons far 
away in the Brillouin zone. It is then convenient to intro- 
duce two continuum fields V&j, one for each pair of nodes, 
expanded around ii^i, ±K2'. 



'i^jaa{q) = Taa{Kj+q). 



(2. 



Here, the wavevectors q are assumed to be small, within a 
circle of radius A around the origin. With this definition, 
the particle/hole transformation is extremely simple. 



(2.9) 



For this reason it is convenient to always define the con- 
tinuum fields 5" around iKj, and account for deviations 
of the node momenta from these values by a particlc/hole 
symmetry-breaking parameter A. 

Once we have restricted attention to the momenta near 
the nodes, it is legitimate to linearize in the quasiparticle 
Hamiltonian. The resulting continuum theory is more 
conveniently written in coordinates perpendicular and 
parallel to the Fermi surface, so we perform the rota- 
tion via X ^ {x — y)/V2 and y ^ {x + y)/-\/2, cor- 
respondingly transforming the momenta and qy (see 
Fig. 2). Linearizing near the nodes, we put eKi+q — vpqx 
where is the Fermi velocity and Ax^+q — ^Qy, where 
A has dimensions of a velocity. An identical lineariza- 
tion is possible around the second pair of nodes, except 
with qx qy Upon Fourier transforming back into real 
space, j{q) — *■ ^j(a?), we arrive at a compact form for 
the Hamiltonian density of the quasiparticle excitations 
in the d-wave superconductor: Tinp = H.^ + TL\ with 



= -^{[vpTHdx + (Ar^ 
+ (1 ^2]x ^ y), 



ydy]^i 



Hx = A^-jT^^-j. (2.11) 



The quasiparticle Hamiltonian takes the form of (four) 
Dirac equations, and can be readily diagonalized giving 
a dispersion relation for the first pair of nodes, 



(2.12) 



and a similar expression with g^; and qy exchanged for 
the second pair. Notice that non-zero A indeed shifts the 
positions of the nodes. 



C. Quantum Fluctuations 

Up to this point, we have taken a spatially constant 
gap function A(x). To disorder the d-wave supercon- 
ductor it is necessary to allow for quantum fluctuations 
of this order parameter. It is tempting to uniformly 
suppress the complex order parameter, and simply put 
A = 0. But doing so recovers the conventional metallic 
state with a Fermi surface. Our task is trickier, since we 
are searching for an intermediate phase, which has strong 
local d— wave pairing (which destroys the Fermi surface) 
but with longer length scale quantum fluctuations de- 
stroying the superconducting phase coherence. Our task 
is similar to the problem of describing the hexatic phase 
in a classical two-dimensional triangular solid,ll3 which is 
intermediate between the crystalline and liquid phases. 
Guided by this example and the principle of pairing be- 
low T* , we want to Gx the magnitude of the complex pair 
field, and introduce fluctuations of its phase. 

Pursuant to this goal, we write 



A 



v^e 



(2.13) 



where is real and (p can be interpreted as the phase of 
the complex superconducting order parameter. The BCS 
gap equation has a degenerate manifold of solutions, for 
arbitrary phase ip. This degeneracy is responsible for the 
Goldstone modes, wherein ip varies slowly in both space 
and time. Our goal is to obtain an effective theory for 
the space and time dependence of if, similar in spirit ta 
the non-linear sigma models "derived" for localization.^ 
Speciflcally, we focus on spatial variations of (p{x) on 
scales longer than A~^. Since (p can vary spatially, care 
is needed in introducing it into the quasiparticle Hamil- 
tonian: 



Ar+idy v^T+e"^/^{idy)e"^/'^. 



(2.14) 



(2.10) 



This symmetric form leads to an hermitian Hamiltonian, 
physical currents, and respects the sy mme tries of the 
problem. A careful derivation of Eq. 2.14 is given in 
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Appendix A. With this prescription, the quasiparticle 
Hamihonian becomes 

+ {1^2;x ^ y). (2.15) 

Since ip can also fluctuate with time, it will convenient to 
consider the time dependence via a Lagrangian formula- 
tion. The Lagrangian density is 



(2.16) 



The appropriate Lagrangian for the phase of the d- 
wave order parameter is simply 



1 



(2.17) 



where the Greek index /i runs over time and two spatial 
coordinates: = 0,1,2 — t,x, y. Here kq is equal to the 
compressibility of the condensate and Kj = —v^kq (for 
j = 1,2 = X, y) with Vc the superfluid sound velocity. We 
expect that the pair compressibility kq is approximately 
one half the electron compressibility of the original elec- 
tron model - in the absence of interactions. If the pairing 
is electronic in origin, one expects that the scale for the 
"charge velocity" Vc is the Fermi velocity. 

As discussed in the introduction, treatment of quan- 
tum phase fluctuations is complicated by the mixing of 
particle and hole variables via the complex gap function. 
To isolate the uncertain charge, we therefore perform a 
change of variables, defining a new set of fermion fields 
ijjj via 



'ipj = exp{-i(pT'' /2)'^j 



(2.18) 



In the superconducting phase, and in the absence of 
quantum fiucutations of the order-paramater phase, one 
can set (p = 0, and these new fermions are simply the 
d-wave quasiparticles. However, when the field if is dy- 
namical and fluctuates strongly this change of variables 
is non-trivial. In particular, the new fermion flelds V' are 
electrically neutral, invariant under a global charge U{1) 
transformation. As we shall see, when the d-wave super- 
conductivity is quantum disordered, these new fields will 
play a fundamental role, describing low energy gapless ex- 
citations, centered at the former nodes. For this reason, 
we refer to these fermions as nodons. For completeness, 
we quote the symmetry properties of the nodon field un- 
der a particle/hole transformation. Since (p —(p, one 
has simply 



(2.19) 



+ il^'2,x^y), (2.20) 
interacting with the phase of the order-parameter: 

Cint ^ df^ipJp. (2.21) 
Here the electrical 3-current is given by 



Jo = l^y^,, 



(2.22) 



(2.23) 



Because the transformation in Eq. 2.18 is local, identical 



expressions hold for these currents in terms of the quasi- 
particle fields, ^E". The form of the particlc/holc asymme- 
try term remains the same in terms of the nodon fields: 



(2.24) 



It is instructive to re-express the components of the 
currents back in terms of the original electron opera- 
tors. One finds 



(2.25) 



(with an implicit spin summation) which corresponds 
physically to the total electron density living at the 
nodes, in units of the Cooper pair charge. Similarly, 



^ ~ 2 3 



-K, 



(2.26) 



corresponds to the current carried by the electrons at 
the nodes. Thus, can be correctly interpreted as the 
quasiparticles three-current. 

To complete the description of a quantum mechani- 
cally fluctuating order parameter phase interacting with 
the gapless fermionic excitations at the nodes, we min- 
imally couple to an external electromagnetic field, A^. 
Since the nodon fermions are neutral, the only cou- 
pling is to the order-parameter phase, via the substitu- 
tion d^Lp — > 9^iy9 — Afj^. For simplicity, here and in the 
rest of the paper, we have set the Cooper pair charge 
2e = 1. The final Lagrangian then takes the form 
C = C^p + + Lint + -Ca, with 



{^^,<^ - Ap)Jp, 



(2.27) 



(2.28) 



The full Lagrangian in the d-wave superconductor, and £0 still given by Eq. 2.20 
C = Cip + Lqp, can be conveniently re-expressed in terms 



of these nodon fields since C 
free nodon piece. 



qp 



C-i, + Cint + C\ with a 



The time component of the electromagnetic field is pro- 
portional to the chemical potential /i, i.e. Aq = 2fj,. For 
electrons at half- filling one has fi — 0. Doping can be 
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achieved by changing fi. Long-ranged Coulomb interac- 
tions could be readily incorporated at this stage by treat- 
ing Aq as a dynamical field and adding a term to the La- 
grangian of the form, Ccoui — (l/2)(9j^o)^- The spatial 
components of the electromagnetic field, Aj, have been 
included to keep track of the current operator. In par- 
ticular, the total electrical 3-current is obtained by dif- 
ferentiating the Lagrangian, i.e. J*°* = dC/dAfj^, which 
gives 



(2.29) 



Here the first terms are the Cooper pair 3-current, and 
the second the quasiparticles current. The equation of 
motion for the phase of the order-parameter, dC/d(p — 0, 
implies the continuity equation d^^Jjf* = 0. 



D. Duality 



(2.33) 



where $ is a (dimensionless) complex field for the 
vortices. The vortex 3-current, following from j'^ = 
dCy/da^, is 



(2.34) 



For small |$| (appropriate close to a second order transi- 
tion) one can expand the potential as, V^{X) = r^X'^ -\- 
u^X"*. The remaining piece of the dual Lagrangian is 



'Ja), (2.35) 



with dual "magnetic" and "electric" fields: b = eijdiQj 
and Cj = v~^{djaQ — doaj). It is straightforward to 
verify that th e du al Lagrangian has the desired prop- 
erty that Eq. 2.32 follows from the equation of motion 
dCo/dafj, = 0. 



To quantum-disorder the d-wave superconductor, one 
must allow for vortices in the pair-field phase, ip. We 
do this using field-theoretic duality, as described, e.g. in 
Ref. To this end we introduce a vortex 3-current, j^, 
which satisfies, 



dyd\ip. 



(2.30) 



In the presence of vortices, is multi- valued, 9^ 1^9 is not 
curl-free, and j"^ is non- vanishing. In the desired dual 
representation, the vortices become the quantized parti- 
cles, rather than the Cooper pairs. However, even in the 
dual representation one still needs to conserve the total 
electrical charge. This can be achieved by expressing the 
total electrical 3-current as a curl, 



J. 



(2.31) 



where we have introduced a "fic titious" d ynamical gauge 
field, a^. Upon combining Eqs. 2.29 - 2.31 , one can elimi- 
nate the pair-field phase, Lp, and relate to the vortices: 



d„[Kx^exapdaap + Ax - Kx^J^] 



(2.32) 



wher e Ju is th e quasiparticle 3-current defined earlier in 
Eqs. ^I2[j2.23| . 

A dual description is obtained by constructing a La- 
grangian, Cd, depending on a^, and j'^, whose equa- 
tion of motion, obtained by differentiating with respect to 
a^, leads to the above equation. To assure that the vor- 
tex 3-current is conserved, it is useful to introduce a com- 
plex field, $, which can be viewed as a vortex destruc- 
tion operator. Since a vortex acquires a 27r phase upon 
encircling a Cooper pair, the vortex field should be min- 
imally coupled to a^. The appropriate dual Lagrangian 
can be conveniently decomposed as £d = £^ + £„ + 
where is given in E(|| |]2^. The vortex piece has the 
Ginzburg-Landau form,@ 



III. NODAL LIQUID PHASE 

In this section we employ the dual representation of 
the d-wave superconductor to analyze the quantum dis- 
ordered phase - a new phase of matter which we refer to as 
a nodal liquid. The dual representation comprises a com- 
plex vortex field, which is minimally coupled to a gauge 
field, as well as a set of neutral nodon fermions. Without 
the nodons and in imaginary time, the dual Lagrangian 
is formally equivalent to a classical three-dimensional su- 
perconductor at finite temperature, coupled to a fiuctu- 
ating electromagnetic field. To disorder the d-wave su- 
perconductor, we must order the dual "superconductor" 
- that is, condense the vortices. The nature of the result- 
ing phase will depend sensitively on doping, since upon 
doping, the dual "superconductor" starts seeing an ap- 
plied "magnetic field" . Below, we first consider the sim- 
pler case of half-filling. We then turn to the doped case, 
where two scenarios are possible depending on jwhether 
the dual "superconductor" is Type I or Type II. □ 



A. Half-filling 

Specialize first to the case of electrons at half-filling, 
with particle-hole symmetry. In the dual representation, 
the "magnetic field" , b, is equal to the deviation of the to- 
tal electron density from half-filling. Thus at half-filling 
(6) — and the dual Ginzburg-Landau theory is in zero 
applied field. The quantum disordered phase corresponds 
to condensing the vortices, setting ($) = $0 7^ 0. In this 
dual Meissner phase, the vortex Lagrangian becomes 



(3.1) 



It is then possible to integrate out the field a^. The full 
Lagrangian in the nodal liquid phase is then 



7 



(3.2) 



t 

k 



(3.4) 



where we have introduced the physical magnetic and elec- At low energies in the nodal liquid phase one can focus 
trie fields: B = e.^d^A^ and = D^Aq - dtAj. The last momenta near the nodes: k = ±K,. The electron 



two terms describe a dielectric, with magnetic perme- 
ability Ho = kq^q and dielectric constant eo = {ijLov'^)~^, 
with the sound velocity entering, rather than the speed 
of light. The external electromagnetic field is coupled to 
the 3-current which can be expressed as a bi-lincar of 
the nodon fermions as, 



eo 



(3.3) 



Notice that this 3-current is automatically conserved: 
0. 

The order {dJ)^ terms which we have not written out 
explicitly are quartic in the fermion fields, and also in- 
volve two derivatives. Since describes Dirac fermions 
in 2 -|- 1 space-time dimensions, these quartic fermion 
terms are highly irrelevant, and rapidly vanish under a 
rescaling transformation. Thus, in the absence of exter- 
nal electromagnetic fields, the description of the nodal 
liquid phase is exceedingly simple. It consists of four 
neutral Dirac fermion fields - two spin polarizations 
(a = I, 2) for each of the two pairs of nodes. 

Despite the free fermion description, the nodal liquid 
phase is highly non-trivial when re-expressed in terms of 
the underlying electron operators. Indeed, the ip fermion 
operators are built from the quasiparticle operators ^E" in 
the d-wave superconductor, but are electrically neutral, 
due to the "gauge transformation" in Eq. 2.18. More- 



over, in the nodal liquid phase, the Cooper pairs are 
not superconducting, but rather in a dielectric Mott- 
insulating phase,ll3 immobilized by their commcnsurabil- 
ity with the underlying crystal lattice. Although the ip 
fermions are electrically neutral, they do carry a new 
conserved "charge". In fact, there are four new con- 
served charges, since the Lagrangian is invariant under 



the global transformations ip 



jaa 



jaa 



for arbi- 
trary constant phases, 9ja, with j = 1,2 and a = 1,2. 
We refer to the -0 fields as nodon operators, their quanta 
as nodons, and the associated conser ved q uantities as 
"nodon charges". As seen from Eq. 2.20, these con- 



served charges are related to the quasiparticle current, 
since Jj = {vp/2)ip'jipj. However, in the nodal liquid 
phase the electrical current operator is not J„ , since 
= dCni/dA^. 



1. Spin response 

The spin response functions in the nodal liquid are 
rather straightforward, since the electron spin operators 
have a simple representation in terms of the nodons. In 
particular, consider the spin operator for small momen- 
tum. 



operators near the nodes can be rewritten in terms of 
the nodon operators, and one finds that back in real 
space the long-wavelength piece of the spin operator, 
S(^) = Y.q cxp(ig • f)Sg, is simply 



S(f) = \ll]]a{x)<Tll^jai^)- 



(3.5) 



Spin correlation and response functions can then be com- 
puted from the free nodon theory. For example, the uni- 
form spin susceptibility is given by 



X 



dE{-df/dE)pn{E), 



(3.6) 



where the nodon density of states is Pn{E) = 
{const)E/vFV&, and f{E) is a Fermi function. One finds 
X ~ T/vpv^. There are also low energy spin excita- 
tions at wavevectors which span between two different 
nodes. The associated spin operators can readily be be 
re-expressed in terms of the nodon fields. For example, 
the staggered magnetization operator, S^, is 



S^f = i [7/;t(T2'crcr2')V't -l-h.c.] 



(3.7) 



Notice that this operator is actually "anomalous" in 
terms of the conserved nodon charge. We will return 
to the effects of finite wavevector magnetic fluctuations 
and ordering in Section IV. 

In addition to carrying spin, the nodons carry energy, 
and so will contribute to the thermal transport. At finite 
temperature, Umklapp scattering processes (or impuri- 
ties) give a finite thermal conductivity; in their absence 
the nodon thermal conductivity is infinite. 



2. Charge response 

The electrical charge properties in the nodal liquid 
phase are, however, somewhat trickier. Imagine chang- 
ing the chemical potential away from /i = 0. In terms of 
the dual vortex "superconductor" this corresponds to ap- 
plying an external "magnetic" field, due to the coupling 



C,. 



-2/i5. The vortices, however, are in the "Meissner" 



phase, and for fJ, < iic the applied field will be screened 
out, maintaining the internal field at 6 = 0. That is, the 
electron density will be pinned at half-filling, until the 
chemical potential passes through the Mott_gap for the 
insulating system of Bosons (Cooper pairs) .Ej 

Despite the presence of a charge gap, there are low en- 
ergy current fluctuations in the nodal liquid. Indeed, in 
this phase the electrical current operator is which is 
bi- linear in terms of the nodon fermions, ^. One might 
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imagine employing this current operator to compute the 
electrical conductivity in the nodal liquid. For this, one 
requires computing a two-point correlator of the elec- 
trical current operator at zero wavevector (say in the 
X— direction) Ix{q = 0) = [e^/ Hov1)diJx{q = 0). But no- 
tice that Jx{q — 0) is proportional to a globally conserved 
nodon charge, since Jx{x) = {vp/2)^\ip^. Thus, when 
the nodon number is conserved one has Ix{q = 0) = 0, 
and the nodons do not contribute to the electrical con- 
ductivity. (There will of course be a response at finite fre- 
quencies in the imaginary part of the conductivity from 
the Mott-insulating phase of the Cooper pairs.) 

When impurities or Umklapp scattering is present, 
however, the nodon number is no longer conserved, and 
the nodons presumably will contribute to the real part of 
the electrical conductivity, at least at finite frequencies. 
Specifically, the umklapp scattering term with momen- 
tum transfer 2?? is given by 



•'Umklapp 



u e 



ABCD 



ll^jAl/jjElpjci'iD + h.C. 



(3.8) 



where the composite index A runs over 1,2,3,4 corre- 
sponding to aa = 11, 12, 21, 22. By power counting, this 
term is irrelevant by one power of frequency. Hence, it 
enters scaling forms in the combination Ueff = uuj. 
According to the Kubo formula. 



cr(a)) ~ 1 (4 (g = 0, uj)lx [q = 0, -uj)) 
^UJ^ {Jx{q^O,uj)Jx{q^O,-Uj)). 



(3.9) 



From scaling, we expect the latter correlation function to 
vary as to. However, as we noted above, it actually van- 
ishes in the absence of Umklapp scattering; therefore it is 
determined by the correction to scaling, which, naively, 
is of the form: 



/ \ 3 2 2 6 



(3.10) 



Finally, it is instructive to consider the behavior of 
the electron Green's function, which can be accessed in 
photo-emission and tunneling experiments. The electron 
operator Ca {x) can be conveniently decomposed in terms 
of the nodon operators by focusing on momenta near the 
nodes. For example, near the node at Kj one can write. 



c„(f) = e^^^-"V^("^)/2V;,ia(f) + 



(3.11) 



In the nodal liquid phase, the electron Green's function, 
G{x,t) = {c^x,t)c{0,0)) factorizes as, 

G(f , t) - e'^^- (e-^'^(-.*)/2e^'^(o" (^t^^ ^)^^.^^ 0)) . 

(3.12) 

Although the nodon correlator is a power law, falling off 
as \x\~'^ and one expects the correlator over exponen- 
tials of the pair field phase to fall off exponentially in the 
nodal liquid, since the Cooper pairs (Bosons) are locked 



in a Mott insulating phase. This indicates a gap in the 
electron spectral function at the Fermi energy, of order 
Ai ~ fic at the nodes. If is small relative to T*, the 
corresponding gap will show strong four-fold anisotropy 
in momentum space, varying from of order T* down to 
of order /ic near the nodon wave- vectors. In the discus- 
sion section, we comment briefly on how such a d-wave 
pseudo-gap feature is likely enhanced when the NL is 
doped in the presence of impurities. 



B. Doping the Nodal Liquid 

We now consider the effects of doping charge into the 
nodal liquid phase. In a grand canonical ensemble this is 
achieved by changing the chemical potential, ^ = Aq/2. 
In the dual Ginzburg-Landau description of the vortices, 
a chemic al pot ential acts as an applied dual field, as seen 



from Eq. 2.35 



smce 



£p = -2fib. 



(3.13) 



The dual magnetic field, b — eijOiUj, is the total electric 
charge in units of 2e. For a hole doping with concentra- 
tion X, one has (6) = x/2a^, with ao the crystal lattice 
constant. Provided the applied dual field, 2/i, is smaller 
than the critical field (2/ic) of the Ginzburg-Landau the- 
ory, the dual superconductor stays in the Meissner phase 
- which is the nodal liquid phase at half-filling. But for 
/i > /ic dual fiux will penetrate the Ginzburg-Landau 
superconductor, which corresponds to doping the nodal 
liquid phase. The form of the dual fiux penetration will 
depend critically on whether the dual Ginzburg-Landau 
theory is Type I or Type II. Within a mean-field treat- 
ment this is determined by the ratio of the dual penetra- 
tion length, A^, to the dual coherence length, (where 
the subscript v denotes vortices). In particular. Type 
II behavior is expected if A„/<^„ > l/\/2, and Type I be- 
havior otherwise. In the Ginzburg-Landau description 
determines the size of a dual fiux tube, which is simply 
a Cooper pair. We thus expect that A„ will be roughly 
equal to the superconducting coherence length, ^, which 
is perhaps 15 — 20A in the cuprates. On the other hand, 
^1, is the size of the "vortex-core" in the dual vortex field, 
and presumably can be no smaller than the microscopic 
crystal lattice spacing, > 3 — 5A. This reasoning sug- 
gest that \v/£,v is probably close to unity, so that either 
Type I or Type II behavior might be possible - and could 
be material dependent. We first consider such Type II 
doping, returning below to the case of a Type I Ginzburg- 
Landau theory. 



1. Type II Behavior 

The phase diagram of a clean thre&dimensional type 
II superconductor is well understood.u Above the lower 
critical field, Hd, fiux tubes penetrate, and form an 
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Abrikosov flux lattice - usually triangular. As the applied 
field increases the flux tubes start overlapping, when their 
separation is closer than the penetration length. Upon 
approaching the upper critical fleld Hc2 their cores start 
overlapping, the Abrikosov flux lattice disappears, and 
the superconductivity is destroyed. Mean field theory 
predicts a second order transition at Hc2, but with ther- 
mal fluctuations one expects this to become weakly flrst 
order .c2l This weak first order transition separates a flux- 
lattice phase from a non-superconducting flux-liquid. 

These results hold equally well for our dual Ginzburg- 
Landau superconductor, except that now the direction 
parallel to the applied field is actually imaginary time. 
Moreover, the Ginzburg-Landau order parameter de- 
scribes quantum vortices, and the penetrating flux tubes 
are Cooper pairs. Upon doping the nodal liquid with 
fi > fj,ci, charge is added to the 2d system, which corre- 
sponds to the penetration of dual magnetic flux. In this 
dual transcription, the resulting Abrikosov flux-lattice 
phase is a Wigner crystal of Cooper pairs, with one 
Cooper pair per real space unit cell of the lattice. We 
denote this paired Wigner crystal phase by PWC. Upon 
further doping, one passes via a weak first order transi- 
tion (at /i = /ic2) into the dual flux-liquid phase. In this 
phase the lattice of Cooper pairs has melted, and they 
are free to condense - this is the d-wave superconductor. 
This latter transition should occur when the spacing be- 
tween dual flux tubes becomes roughly comparable to the 
coherence length, ^y. Experimentally, superconductivity 
typically sets in for x = 0.1, which corresponds to one 
Cooper pair for every 20 or so Cu atoms, and a mean 
pair separation of 4-5ao. This again suggests that Xy/Cv 
is probably of order one. 

In the Cooper pair Wigner crystal phase, translational 
symmetry is spontaneously broken. However, in a real 
material the Wigner crystal will have a preferred loca- 
tion, determined by impurities and perhaps crystal fields, 
which will tend to pin and immobilize the Wigner crys- 
tal. The resulting phase should be an electrical insu- 
lator. Moreover, in two-dimensions even weak impuri- 
ties will smooth the weakly first order transition between 
the Wigner crystal and superconducting phases. (In the 
absence of impurities, long-ranged Coulomb interactions 
preclude phase separation, so a mixed phase would result 
- see Type I behavior below.) 

A striking and unusual feature of the PWC phase, is 
that it co-exists with the nodal liquid, as we now ar- 
gue. With a weak (commensurate) pinning potential 
present, the Wigner crystal phase is a dielectric. The 
charge response is thus essentially the same as that of 
the undoped phase at half-filling, except with a modi- 
fied dielectric constant and magnetic permeability. Thus 
even with doping, it is possible to "integrate out" the 
charge fiuctuations described by the fields a^, and ar- 



rive at the nodal liquid Lagrangian Eq. 3.2, except with 



different values of eo and /io- The only complication is 
that there will be a background frozen in charge den- 
sity, from the Wigner crystal of Cooper pairs, so that 



{b{x)) = n{x) = no + 5n(x). Here, the mean pair density 
is simply uq — x/2aQ for doping x, and Sn(x) has the pe- 
riodicity of the Wigner crystal. This "ba ckground" field 
couples to the nodons, and from Eq. ^.35 leads to a term 
of the form, 

A = {b{x))i^^'j„ = -Ln{S-)^y^Pj. (3.14) 

What is the effect of this term on the nodons? Consider 
first the spatially constant piece, proportional to no. This 
term can be absorbed into which contains a term of 
the form, ip^VFqjT^ipj, and leads to a momentum space 
shift of the nodes, with qj qj + (no/2wi?Ko). Since 
the pair compressibility kq ^ dno/dfi, the shift satisfies 
vpSq Sfi, as expected from the change of the area en- 
closed by the Fermi surface upon doping. 

The spatially varying background density, 6n(x), 
which has the periodicity of the underlying Wigner crys- 
tal, causes a mixing between nodon states at momen- 
tum differing by Wigner crystal reciprocal lattice vectors. 
When the reciprocal lattice vectors are larger than the 
momentum cutoff. A, of the nodons, the dn{x) term can- 
not scatter within the low energy nodon theory, and can 
be dropped. At lower pair densities, it becomes necessary 
over a range of length scales to retain the new periodicity 
by working with nodon Bloch states, rather than plane 
waves. (Since only linear derivatives enter into the nodon 
Lagrangian, Bloch wavefunctions in the lowest band can 
be readily constructed.) In any event, at length scales 
larger than the Wigner crystal lattice spacing, the form 
of the effective theory of the nodons is identical to that 
at half- filling. 

We thereby arrive at a description of a rather remark- 
able new phase of matter. A Paired Wigner Crystal 
(PWC) of doped Cooper pairs co-exists with neutral 
gapless fermionic excitations - the nodons. In this co- 
existing phase, which we denote as PWC/NL, low en- 
ergy spin and thermal properties will be dominated by 
the nodons. The behavior will be qualitatively similar to 
that in the undoped nodal liquid phase. We propose that 
this PWC/NL phase is present in the pseudo-gap region 
of the high Tc cuprates. 



2. Type I Behavior 

In a classical Type I superconductor, thg applied field 
is expelled until the critical He is exceeded.B At this point 
there is a first order phase transition from the Meissner 
phase with all the flux expelled, to a normal metal phase 
in which (essentially) all the fleld penetrates. If a thin 
film type I superconductor is placed in a perpendicular 
field, screening currents are unable to expel all the flux, 
and a "mixed" or "intermediate" state occurs. In this 
mixed phase, regions of superconductivity co-exist with 
normal metallic regions. In some cases the superconduct- 
ing regions form stripes, but generally the lowest energy 
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configurations are determined in large part by material 
imperfections, and tend to be "history" dependent. 

If our dual Ginzburg-Landau theory describing quan- 
tized vortices is of type I, then similar properties are ex- 
pected. Specifically, as the chemical potential increases, 
the dual field - which is the Cooper pair density - re- 
mains at zero until a critical chemical potential /ic is 
reached. At this point there is a first order phase tran- 
sition, between the nodal liquid phase at half-filling, and 
a d-wave superconductor at finite doping, Xc- At fixed 
doping X < Xc, phase separation is impeded by long- 
ranged Coulomb interactions between the Cooper pairs. 
The system will break apart into co-existing "micro- 
phases" of nodal liquid and d-wave superconductivity. 
The configuration of the co-existing "micro-phases" will 
be determined by a complicated competition between the 
Coulomb energy and the (positive) energy of the domain 
walls. In practice, impurities will also probably play a 
very important role. This doping scenacio is similar to 
that envisaged by Emery and Kivelson,tj who have ex- 
tensively discussed the possibility of phase separation as a 
mechanism for high superconductivity. Unfortunately, 
with the transition being strongly first order in this case, 
the associated physics is rather non-universal. 



IV. ANTIFERROMAGNETISM IN THE NODAL 
LIQUID 

A. Effective Action for Antiferromagnetism 

We now turn to the low-doping region of the phase 
diagram of Fig. 1. Retaining the Nodal Liquid as the 
underlying description of the low-energy fermionic de- 
grees of freedom, we consider antiferromagnctic order- 
ing. In principle, this can arise in two ways. Antiferro- 
magnetism could stem from interactions between nodons, 
i.e. physics below the scale T* . This can be modeled 
in principle by including simple inter-nodon interactions. 
The experimental coincidence of T* with the Neel tem- 
perature and magnon bandwidth J suggests that such 
a separation of scales is not valid. Instead, antiferro- 
magnctic correlations may exist already at (high) en- 
ergies comparable to T* . This sort of local AF am- 
plitude could be captured by decoupling spin-spin (or, 
e.g. on-site Hubbard) interactions in a microscopic model 
with a Hubbard-Stratonovich transformation. Such a 
decoupling introduces a conjugate field M, which inter- 
acts with the electron operators via a term of the form 
Hm = Z]sM(f)-c^cra^c^. Integrating out high-energy 
degrees of freedom generates an effective action for M. 
We expect dominant ordering tendencies at momentum 
7? = (7r,7r), and so decompose 



M Mo -I- exp(i7f • x)N, 



(4.1) 



uniform and staggered magnetization. Focusing on the 
Neel ordering, we imagine integrating out Mq to obtain 
the Lagrangian 
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C = -i^^|9,,N|^ - Vn{\N\) + .gN • S^, 



(4.2) 



where Kq = if, Ki = K2 = —v^K, with Vs the spin- 
wave velocity in the AF. Here S^f is the spin operator 
at momentum tt, expressed as a bi-linear in terms of the 
electron operators as in Section III. The staggered mag- 
netization operator can be readily re-expressed in terms 
of the nodons as, 



(4.3) 



where N and Mq are slowly- varying. The fields Mq and 
N have the physical interpretation of the coarse-grained 



given earlier in Eq. 3.7. Near any phase transitions, 
and for most phenomenological purposes, it is sufficient 
to take a simple form for the potential: VVdA^I) = 
?'Af|A^P + un\N\*. The parameter r^r controls the pres- 
ence or absence of AF order. In mean-field theory, and 
neglecting for the moment the nodon coupling g, the 
ground state passes from long-range to short-range AF 
order as tat is tuned from negative to positive. A precise 
determination of r^y, uat, and g in terms of i, t' , U, x, 
etc. is the province of a more mic rosc opic theory. 

We note that, in principle, Eq. 4.2 allows for the pos- 
sibility of incommensurate spin-density-wave ordering at 
wavevectors other than (tt, tt), which would correspond 
to a state which spontaneously develops a spatially peri- 
odic expectation value for N. We find such ordering un- 
likely within this model, however, and have therefore not 
included terms responsible for locking in possible higher- 
order commensurate magnetic wavevectors. Since incom- 
mensurate order seems not to be realized experimentally 
at low doping, we hope this omission is unimportant. 



B. Magnetism and Phases at Half-Filling 

Once we have coupled in the Neel order parameter 
field, we can describe magnetic phases, in addition to 
the nodal liquid and d-wave superconducting phases of 
earlier sections. Here we first focus on the situation at 
half-filling, where our effective field theory already de- 
scibes a number of magnetic and non-magnetic ground 
states. It will be useful to further specialize initially to 
models with particle-hole symmetry, returning later to 
the half-filled but particle-hole asymmetric case below. 



1. Particle-hole symmetric case 

The full effective Lagrangian has two order parameter 
fields, the Neel order parameter, N, and the vortex com- 
plex field, <&, which is minimally coupled to a gauge field, 
a^. The Neel order parameter is directly coupled to the 
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nodons, whereas the vortex field only sees the nodons in- 
directly via the gauge field. Ordering of the two fields is 
determined by the coefficients of the quadratic terms in 
the Lagrangian, namely rjy and r$. It will be convenient 
to plot the phase diagram at half- filling in the r^r — r$ 
plane. The phase diagram with particle-hole symmetry 
(A = 0) is shown in Fig. 3a. Here we briefly discuss each 
of the four phases. 
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Fig. 3: Phase diagrams at half-filling for the particle- 
hole (a) symmetric and (b) asymmetric cases. In (b), 
both horizontal phase boundaries shift downward with 
increasing particle/hole asymmetry parameter A, dimin- 
ishing the domain of the AF phase. 

Consider first large and positive, so that Neel order 
is not present. For r$ negative, the vortices will con- 
dense (the Meissner phase in the dual Ginzburg-Landau 
theory) leaving the nodons as the only low energy gap- 
less excitations. This is the nodal liquid phase. As r$ 
changes sign, the vortices will disorder, entering the "nor- 
mal" phase of the dual Ginzburg-Landau theory. This 
phase corresponds to the d-wave superconductor, shown 
in the upper right quadrant of Fig. 3a. The d-wave super- 
conducting phase can be obtained in a microscopic lattice 
model even at half-filling, hv appropriately choosing the 
electron interaction terms.E2l 

When rpf is large and negative, the model magneti- 
cally orders into the antiferromagnctic Neel phase. With 
particle-hole symmetry the antipodal nodes are separated 
by the Neel ordering wavevector, (tt, tt), so the nodons are 
"nested" . This opens a gap in the nodon spectrum, as 
we now demonstrate. With either sign of r$ the gauge 
field can be integrated out, generating irrelevant four 
fermion nodon interaction terms, which can be dropped. 
Then, upon putting (N) = iVpy into the effective La- 
grangian, we arrive at the following quadratic Lagrangian 
for the nodon fields, 



(4.4) 



This model can be readily diagonalized with an appro- 
priate Bogoliubov transformation, giving energy eigen- 
values, 

EM = ±^J{vFq:,Y + {v^qyY + [gNo^, (4.5) 



in the j = 1 sector, and an identical form with and 
qy interchanged for the other pair of nodes (j = 2). 
In all nodon sectors there is a non-zero gap, equal to 
gNo. In the lower left quadrant of Fig. 3a, with r$ nega- 
tive, this corresponds to the usual Neel antiferromagnet. 
With the nodons gapped out, the only low energy exci- 
tations are the spin waves of the antiferromagnet. For 
r$ positive, in the lower right quadrant, antiferromag- 
netism co-exists with d-wave superconductivity. In the 
d-wave superconductor, the nodons become equivalent to 
the d-wave quasiparticles, so that the d-wave state in this 
quadrant is rather unusual. In particular, it is a d-wave 
superconductor with a full single particle gap, and an ab- 
sence of nodal quasiparticles. We have denoted this with 
a subscript g ~ for gap ~ in the phase of the lower right 
quadrant. 

Before turning to the effects of particle-hole asymme- 
tries, it is interesting to briefiy discuss the nature of the 
phase transitions between the four phases in Fig. 3a. 
Consider first the vertical phase boundary, separating 
the superconducting from non-superconducting phases. 
For rjv negative there are no gapless nodons, and the 
magnet is ordered in both phases. At the transition, 
r$ = 0, where superconductivity develops, we can em- 
ploy the dual Ginzburg-Landau theory. Equivalently, we 
can return to the original representation (before duality) 
in terms of the pair field phase, ip, which_(in the ab- 
sence of long-ranged Coulomb interactionfO) is simply 
the classical three-dimensional XY model. The result- 
ing transition is in the classical 3d-XY universality class. 
For rpj positive, we have to worry about the presence of 
gapless nodons, which might effect the nature of the su- 
perconducting transition at r$ = 0. However, since the 
dual Ginzburg-Landau field $ is only indirectly coupled 
to the nodons, via the gauge field a^, we expect that the 
transition will still be in the 3d-XY universality class. In- 
deed, power counting in d spatial dimensions about the 
Gaussian fixed point, reveals that the coupling term (of 
the form datp''il;) has scaling dimension {3d+l)/2. Being 
greater than d + I for all d > 1, and much greater near 
the upper-critical dimension (due = 3), one expects this 
coupling to be strongly irrelevant in the two-dimensional 
case of interest. 

Consider next the phase transition into the antiferro- 
magnetic state, upon crossing the horizontal axis. For 
either sign of r$, since the vortices and gauge field 
can be integrated out generating irrelevant four-nodon 
interactions, the transition is described by the Neel order- 
ing field coupled to the nodon bi-linear, with Lagrangian 
C ~ £0 + Cn- This is an interesting field theory de- 
scribing Dirac fcrmions anomalously coupled to a fiuc- 
tuating 0(3) field. Power counting about the Gaussian 
fixed point (free nodons and rN — un = 0) reveals that 
the coupling term is relevant in two dimensions, but if 
the model is suitably generalized to higher dimensions 
becomes marginal in d = 3. This suggests an attack near 
four space-time dimensions, working perturbatively for 
small e = 4 — (d + 1). A complication is that there are 
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in fact three independent velocities, two for the nodons 
(vf,v^) and one for the Neel spin field {vs). In Appendix 
B we first consider the special case in which all three ve- 
locities are set equal. This model is then a bona fide rel- 
ativistic field theory, with Lorcntz invariance. A leading 
order perturbative renormalization group in e reveals the 
existence of a new non-trivial fixed point. The associated 
critical properties and exponents are briefly discussed in 
Appendix B. Remarkably, as we also show in Appendix 
B, this relativistic fixed point is in fact linearly stable to 
small deviations in the three velocities. A microscopic 
model with different velocities scales into a relativistic 
form near criticality. As further discussed in Appendj^ 
B, this model is closely related to a remarkable modeled 
which has has a non-trivial fixed point in D = 2-4- 1 with 
an exact (TV =1) space-time supersymmetry.EJ Super- 
symmetry is powerful enough-tp determine several exact 
critical exponents in _D — 3p3 which serve as a useful 
check on our e expansion results. 

It is finally worth mentioning the effects of the coupling 
to the nodons, on the location of the antiferromagnetic 
ordering transition. With positive one can safely in- 
tegrate out the nodons. This leads to a supression of tat, 
with a "renormalized value", fjv, given to leading order 

by, 



(4.6) 



Here A is the high momentum cutoff on the nodons, and c 
a cut-off dependent constant of order one. The coupling 
to the nodons thus tends to enhance the magnetic or- 
dering. Within mean-field theory, the antiferromagnetic 
ordering transition will take place when the "renormal- 
ized" coefficient fjv = 0, as indicated in Fig. 3a. 



2. Effects of Particle-hole asymmetry 

Staying at half-filling, we next consider the effects 
of particle-hole asymmetries, generated for example by 
a second neighbor electron hopping term in a micro- 
scopic square lattice model. As discussed in Section II, a 
particle-hole asymmetry generates an additional term in 
the nodon Lagrangian of the form. 



(4.7) 



In the absence of coupling to the Neel field, this sim- 
ply causes a momentum space shift in the positions of 
the nodes, by an amount 6q = X/vp- As we shall see, 
this has a profound effect when the nodons are coupled 
to the Neel order parameter, since such a shift destroys 
the "nesting" of the nodons. Indeed, this leads to two 
additional phases at half-filling, as shown in the phase 
diagram for A ^ in Fig. 3b. 

Since both additional phases are antiferromagnetically 
ordered, we once again integrate out $ and (valid 
for non-zero r$), and put (N) — Noy into the effective 



Lagrangian to arrive at a quadratic nodon Lagrangian 
of the form Cnodon + with Cnodon given in Eq. 4.4. 



Once again a Bogoliubov transformation diagonalizes the 
quadratic form, and with non-zero A the energy eigenval- 
ues (in the j = 1 nodon sector) satisfy, E'^—2AE^-\-B = 
with, 



and 



A = [gNof + [vpq.f + (v^qyf + A^, (4.8) 



B =[A-2X^] + {2Xv^qyf. (4.9) 



If gNo > A, there is no solution at = 0, so that there is 
a gap in the nodon spectrum. The resulting phases - an- 
tiferromagnet for r$ < and co-existing gapped d-wave 
superconductivity with antiferomagnetism for r$ > - 
occur for large negative rjv, where Nq is large (given by 
Nq = ■\/ (— rAr)/2ujv within mean-field theory). But for 
smaller |rjv|, when gNo < A, zero energy solutions do 
exist, and there are gapless nodon states present! This 
leads to the two new phases present in Fig. 3b. 

Specifically, for r$ negative, the new phase exhibits 
gapless nodon excitations co-existing with long-range an- 
tiferromagnetic order. The nodons are now incommensu- 
rate with the magnetic order (at tt ) , since the zero energy 
nodon state for A > gNQ occurs at a shifted wavevector: 



l. = Vp'^X^~igNor^O, 



(4.10) 



and qy = 0. This interesting new phase, which we de- 
note as AF/NL, exhibits gapless incommensurate mag- 
netic fiuctuations co-existing with the tt magnons of the 
Neel state. 

For positive r$, where d-wave superconductivity is 
present, the gapless nodon excitations are simply d-wave 
quasiparticles, incommensurate with the Neel order. In 
the new phase, the incommensurate gapless d-wave su- 
perconductor co-exists with antiferromagnetic order, as 
depicted in Fig. 3b 

It is worth emphasizing that with particle-hole asym- 
metry present (ie. non-zero A) the region of antiferromag- 
netic order in the phase diagram is diminished, relative 
to the case with A = 0, as depicted in Fig. 3b. This 
occurs because incommensurate nodons are less effective 
at renormalizing r^. Specifically, upon integrating out 
the nodons with non-zero A one finds the same form as 
in Eq. 4.6, but with c cF{X/vfA), where F{X) is a 
monotonically decreasing function of X with F{0) = 1. 
This leads to a downward shift in the horizontal phase 
boundary in Fig. 3b by an amount. 



6fN - [1 - FiX/vFA)]{cA/v^)g^ 



(4.11) 



Physically, particle-hole asymmetries, such as a second 
neighbor hopping term, tend to frustrate and weaken the 
antifcrromagnetism at tt. With non-zero A the nodes are 
shifted off commensurablilty, and the magnetism is in- 
deed weakened. This effect leads to a natural mechanism 
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for the destruction of antiferromagnetism upon doping, 
as we describe in the next section. 

Finally, we briefly discuss the nature of the phase tran- 
sitions between the six phases present at half-filling with 
particle-hole asymmetry. Arguments as in the previous 
subsection, suggest that the vertical phase boundary sep- 
arating the three superconducting from the three non- 
superconducting phases (at r$ — 0) should, as before, 
be in the universality class of the classical Sd-XY model. 
Since the nodons are incommensurate at the upper hor- 
izontal phase boundary where antiferromagnetism first 
appears, they will decouple from the critical magnetic 
fluctutations. Both of these two transitions (for posi- 
tive and negative r$) should thus be in the universality 
class of the classical 3d Heisenburg model. At the lower 
horizontal phase boundaries, the gapless nodons disap- 
pear. The critical properties are correctly described by 
the quadratic nodon Lagrangian, considered above. In 
particular, for A > gNg, one can linearize for small mo- 
mentum around the shifted zero energy nodes by putting 
Sqx = Qx- Vp^^X^ - {gNoY, which gives, 




v&qy 



(4.12) 



As the transition is approached, the velocity along the 
x-direction - i.e. perpendicular to the Fermi surface - 
vanishes and the nodons become quasi- ID. 



two alternative phases in undoped cuprates. Since by 
assumption we are varying only fi, and /i couples indi- 
rectly to m agnetism via the corresponding dual "flux" 
in Eq. 2.35 , the AF order and nodons are effected only 
once charge is added to the system. What is the effect 



of such charge doping on the AF? From Eqs. 2.24 , 2.35 
we see that the added charges act simply to increase the 
efFective particle/hole asymmetry of the nodons. In par- 
ticular, 



^eff - A 



2 • 



(4.13) 



Eq. 4.13 is an extremely useful result. Using it, we can 
simply trade the doping x for an effective particle/hole 
asymmetry to determine the fate of the system from the 
phase diagrams at half-filling. Fig. 3. 
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C. Doping the Antiferromagnet 

We are now in a position to extend our understand- 
ing of doping to include AF order at half-filling, in ac- 
cord with experimental observations. In all high-Tc ma- 
terials, doping is actually achieved by chemical substitu- 
tion/depletion of atoms between the Cu02 layers. While 
it is generally believed that this process transfers charge 
to the Cu02 planes, this charge transfer is not necessarily 
proportional to the chemical doping, defined as the frac- 
tion of atoms substituted or removed. To simplify this 
discussion, however, we shall assume that chemical dop- 
ing indeed corresponds to hole doping, and consider the 
phase diagram as a function of hole concentration x. A 
further assumption concerns the degree of particle/hole 
asymmetry. Since the composition is changing with dop- 
ing, the parameter A should in general be a function of 
X (in fact, we expect the asymmetry to increase with x). 
To simplify the discussion, we shall further assume that 
any explicit dependence of A on a; is weak, and therefore 
treat the effects of doping solely through the chemical 
potential /i. 

At half-filling, the system sustains long-range AF or- 
der. In general, we expect a non-zero particle/hole asym- 
metry, so that this could correspond to either the AF or 
AF/NL state, the latter occuring if A is sufficiently large. 
We do not believe current experiments distinguish the 



Fig. 4: Possible phase diagrams as a function of dop- 
ing for the type II (a and b) and type I (c) scenarios. 
We tentatively identify the PWC/NL phase, which ex- 
hibits neither magnetism or superconductivity, with the 
pseudogap state of the underdoped high-Tc materials. 

Consider then first the type II doping scenario. The 
charge behavior is similar to that obtained when dop- 
ing the NL, Sec. [II. Upon increasing fj, from zero, the 
dual "flux" is first expelled from the sample, and the 
system remains undoped. Charge first enters above the 
dual "lower critical field" , /i > ^ci , forming a Paired 
Wigner Crystal (PWC) with density x{^) due to long- 
range Coulomb interactions. Since x is small at this 
point, the crystal coexists with the AF, so the actual 
phase for small x is an AF/PWC. As x increases, so does 
Xeff, unbinding the nodons into the AF/PWC/NL. This 
can be understood from the evolution of the phase dia- 
grams at half-filling as a function of A, as shown in Fig. 3. 
As X increases further, the NL and AF become increas- 
ingly incommensurate, and the energy gain from their 
coupling is eventually reduced sufficiently to destroy the 
AF order in a transition to a PWC/NL phase - again 
a feature of the phase diagrams at half-filling. Finally, 
when X > Xc2, the upper critical field is reached and 
the crystal melts into the dSC phase. This progression 
is shown schematically in Fig. 4a. An alternate type II 
doping scenario, shown in Fig. 4b, is that the system is an 
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AF /NL at half- filling, in which case the phase diagram is 
unchanged except for the absence of the AF/PWC phase. 

Another possibility is type I doping. Because this 
involves a strong first-order transition in the absence 
of Coulomb interactions, the mixed (micro-phase sepa- 
rated) state could occur as a coexistence between a num- 
ber of different phases. The simplest phase diagram in- 
cludes only coexistence between the AF and pure dSC 
as in Fig. 4c. Because the physics of the mixed state is 
highly non-universal, we do not discuss it further here. 



V. DISCUSSION 

The main result of this paper is the Lagrangian, 
Eqs. |l.l[ - |l.2| , which describes the Nodal Liquid phase, 
its interaction with external electromagnetic fields, and 
transitions between it and the antiferromagnetic (AF) 
and superconducting (dSC) phases. This Lagrangian fol- 
lows directly from disordering the d-wave superconduc- 
tor. It implies that, in the underdoped region, low-energy 
fcrmionic degrees of freedom arc located solely at four 
isolated (Dirac) points in the Brillouin zon&.^a hypoihe.- 
sis which is stronglyj^unported by ARPESQfl, NMRHO, 
optical conductivitytllE3, an d other experimental. The 
main consequences of Eqs. |l.l[ - |l.2| are (1) the predic- 
tion of a new zero-temperature phase, the Nodal Liq- 
uid, which interpolates between the AF and dSC phases; 
(2) a quantitative description of charge and spin dynam- 
ics in this phase; (3) specific predictions for the critical 
behavior at the AF and dSC ordering transitions; and, 
above all, (4) a coherent weak-coupling framework - with 
the Nodal Liquid as its foundation - for understanding 
the underdoped side of the high-Tc phase diagram. We 
tentatively identify the coexisting Paired Wigner Crys- 
tal/Nodal Liquid (PWC/NL) phase as the pseudo-gap 
state in a hypothetical disorder-free underdoped cuprate. 
In the real materials, however, disorder will play a role, 
as we briefly discuss below. 

Our description of this part of the phase diagram en- 
joys kinship with several other approaches. The Nodal 
Liquid phase is reminiscent of the prfiux stateEj and the 
SU{2) MFT staggered-fiux state,c3 and is the d-wave 
analog of the-ahort-rangeEj resonating valence bond spin- 
liquid state.ta These states also involve neutral Dirac 
fermions interacting with a gauge field, but the cou- 
pling to electromagnetic fields, the coupling between the 
fermions and the gauge fields, the bosonic charged de- 
grees of freedom, and the evolution with doping aiM- 
rather different in the Nodal Liquid. Several author£ 
have conjectured that the lightly-doped 3-leg Hubbard 
ladder might serve as a paradijgm for the underdoped 
cuprates and Furukawa and RiceEj have tried to substan- 
tiate these claims with weak-coupling RG calculations on 
partially nested Fermi liquids. The Nodal Liquid con- 
cretely realizes the attractive features of this proposal. 
The basic idea of bringing AF and dSC under the same 



rubric, which is expressad in (^]^), is the central theme of 
Zhang's 50(5) theoryEj However, there is a direct tran- 
sition from AF to dSC in the 5*0(5) theory, whereas the 
Nodal Liquid intervenes in our theory. There is a further 
important distinction, namely, that our theory focuses 
on the zero-temperature quantum phase transitions of 
the high-Tc materials. This is one reason why our theory 
accords primary importance to the low-energy fermionic 
degrees of freedom. Finally, our prediction of phase sepa- 
ration at the dSC transition in the type I scenario as well 
as our intespretation of T* echoes the ideas of Emery and 
Kivelson.ll3 

There are a number of important issues which we have 
not addressed in this paper. By restricting our atten- 
tion to the region underneath the dashed line in Fig. 1, 
we have skirted one of the most controversial questions 
in this field: what mechanism drives pair formation at 
this scale? Presumably this physics must be understood 
for progress to be made on the part of the phase dia- 
gram above and to the right of the dashed line. This 
would require an investigation complementary (but per- 
haps orthogonal in spirit) to ours. Also, our discussion 
of transport was necessarily incomplete because finite- 
temperature transport can be particularly subtle (for re- 
cent examples of this, see Ref. 33) and the transport prop- 



erties of the Nodal Liquid deserve a thorough exposition 
of their own, which we defer. Moreover, a pure sample 
would melt at a finite temperature phase transition, al- 
though this transition would be rounded by impurities. 
The melting temperature is expected to vanish upon ap- 
proaching either zero doping or the PWC/NL to dSC 
quantum phase transition, and therefore has maximum 
at some intermediate x in the underdoped regime. 

The effects of disorder are also quite subtle, and war- 
rant a full and separate treatment. Nevertheless, a few 
comments are germane to this discussion. The first, and 
most basic, is that disorder plays a significant role in the 
physics of the cuprates. Even in YBa2Cu:i07-s, which 
is believed to be cleaner than, say, La2-xSrj.CuOi, 
doping cannot help butpLptroduce disorder. According 
to standard arguments)^ first-order phase transitions 
will be driven second-order by arbitrarily weak disor- 
der in two dimensions. In particular, we expect the 
PWC/NL^dSC transition to be second order with im- 
purities present. Moreover, based on the irrelevance (in 
the technical sense) of the coupling between the nodons 
and the superconducting phase in the clean case, we 
suspect this transition may be in the same universality 
class as the superconductor-insulator transition. This 
could explain the expjeriments of Fukuzumi, et al. on 
Y Ba2Cu3^yZnyO'/-s^ The disorder will also have an 
effect on the phases themselves. For instance, transport 
and spin dynamics in the Nodal Liquid will be influ- 
enced by disorder. Finally, we note that disorder will 
transform the PWC into a Bose Glass (BG). One con- 
sequence would be power-law suppression (rather than 
a hard gap) of the low frequency electron spectral func- 
tion at the nodes in the BG/NL phase because the BG is 



15 



compressible. t3 Having adopted a panoramic view in the 
preceeding section, we can afford, in closing, to narrow 
our focus to the BG/NL (and the PWC/NL from which 
it descends) because it is our candidate for the T — 
pseudo-gap phase: a phase without a Fermi surface or 
long-range order but possessing low-energy fermionic ex- 
citations centered about four points in the Brillouin zone. 
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APPENDIX A: MICROSCOPIC APPROACH 

In this appendix, we describe techniques to derive the 
effective field theory from some specific microscopic mod- 
els. As our purpose is primarily phenomenological, we 
will consider one of the simplest models which develops 
antiferromagnetism and d-wave superconductivity. This 
is a square lattice extended Hubbard model with nearest- 
neighbor hopping t, on-site electron-electron repulsion U 
and nearest neighbor attraction V (we emphasize that 
this model is not realistic, but is chosen for illustrative 
purposes), i.e. 

H[c\c] = -tY, [4(2^^(2?') + h.c.] +M^n(f) 

{xx') X 

+U [n{x)f -VY^ n(f)n(x'). (Al) 

{xx') 

Here n{x) — cl^{x)c^{x), and jj, is the microscopic chem- 
ical potential, and we have neglected to include terms 
(e.g. a second-neighbor hopping t') which break parti- 
cle/hole symmetry at half-filling. Also, for simplicity we 
measure distances in units of the lattice spacing. The 
tendency toward AF and dSC states can brought out by 
using the identities [n(x)]^ — — 4/3[S(a;)]^ + 2ji{x) and 
n{x)n{x') — cl^{x)c'f^{x')Cf^{x')c^{x), for x ^ x' . The 
Hamiltonian can be rewritten as 

B\c\c\^~tY^ [4(f)c„(f' ) + h.c.] 

{xx') X 

-^[/X] |S(f)P " y X] ci{x)cl{x')cp{x')cjx), (A2) 

X {xx') 

where we defined a shifted chemical potential /i, and have 
neglected an unimportant constant. As usual, the lattice 
spin operator is defined by S(x) = icJ,(a;)(Ta^c^(x). The 



angular brackets {xx') under the two sums indicate sums 
over all nearest-neighbor pairs of sites. 

To derive an effective field theory, it is convenient to 
use an imaginary time path integral formulation. In this 
case one studies the partition function Z = Tre^^^^, 
where T is the temperature. It can be represented using 
Grassman coherent states as 



[dc] [dc\e 



where the Euclidean action is 



S 



Jdr i^^Ca{x)drCaix) + H[c,c] 



(A3) 



(A4) 



We consider here only T = 0, for which the r integration 
domain is infinite. The possibility of dSC and AF phases 
can be entertained by decoupling the above action using 
Hubbard- Stratonovich transformations. One finds that 



Z = / [dc][dc][dM][dA][dA]e' 



-Si 



(A5) 



with Si = J dT[J2s (^a{x)drCa{x) + H^ff]. The effective 
Hamiltonian can be decomposed into i?e// = Hqp+HM + 
Ho + H,nt, and 



H^, with Hqp 



Ha 

Hint 



-t [4(2^)0, (f') + h-c.] +/i^n(f) 

{xx') X 

-X]M(f).S(f) 

X 

" Yl ['^ss'C0ix')cc{x) + Ag§,cc{x)c/3{a 

{xx') 



H 



M 



X 

H. = ^Y^'^I'- 



(A6) 

(A7) 
(A8) 
(A9) 



{xx') 



Eqs. A6-A9 form a basis for studying the original ex- 
tend ed H ubbard model. Following the philosophy of 
Sec. HA , we imagine integrating out high-frequency 



modes in the functional integral to arrive at an effec- 
tive field theory for the low-lying degrees of freedom. In 
the process, one will generate dynamics for the order pa- 
rameter A and the magnetization M. For the most part, 
symmetry considerations require the corresponding La- 
grangians to take the forms given in Sec. H A and Sec. |rv| , 
so we choose not to complicate the presentation by ex- 
plicitly performing these integrations (e.g. diagrammat- 
ically). 

One subtle point in the analysis of Sec. 
ever, does warrant a more careful treatment. This is the 
coupling of the nodons to the superconducting phase- 
gradient, from which follows the expressions for the 



HA, how- 



quasiparticle current, Eqs. 2.22 



2.23. In Sec. HA, we de- 



rive these using the "symmetric" prescription of Eq. 2.14 
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We now show that the currents are indeed obtained cor- 
rectly using this prescription. 

We first specialize to the case of singlet pairing, A"-, = 
A{x, x'){SaiS[}i —SaiSpi). Since A lives on the bonds, it 
is convenient to associate two such fields with each site 
in the square lattice, i.e. 



As expected, the 9 field is massive, and can be integrated 
out. Equivalently, one minimizes Ti.2 with respect to 9. 
This process restores isotropy and gives 

7, 



Ho 



(A15) 



Ai(f) = A{x,x + ii), 
A2(f) = A{x,x + e2), 



(AlO) 
(All) 



This is just the spatial component of the superfluid La- 
grangian, Eq. p. 17 , with 7 — v^kq. The corresponding 



where ei, 62 are unit vectors along the a and b axes of the 
square lattice, respectively. Note that at this point, we 
have specified no particular relation between Ai and A2, 
so that the model has the potential both for d-wave and s- 
wave pairing. To distinguish them, we must consider the 
form of the effective action for A, A generated upon in- 
tegrating out the high-energy modes. By symmetry, the 
simplest local allowed additional term on the lattice is 
a sum of /7(l)-invariant two-bond products around each 
lattice site, which can be written as S2 = Jdr H2, with 

H2 = J^|Ai(f)A2(x) + A2(f)Ai(x-ei) 



+Ai{x~ ei)A2{x- 62) + A2(x - e2)Ai(f) + c.c. 



(A12) 

Of course, the actual quadratic AA interaction terms 
generated from the high-energy degrees of freedom will 
be much more complex. However, since the general form 
of the long- wavelength effective action is dictated by 
symmetry, this example su ffices for illustrative purposes. 
Writing A^ = Aoe*'^^ Eq. ^l| becomes 



time component can be obtained similarly. 

Finally, consider the coupling of of the phase (p to the 
fermionic quasi-particle operators. To study this, we 
take for simplicit y /i = M = 0. Using the definitions 
in Eqs. AlO- All, the coupling term in Eq. A7 can be 
rewritten 



Hi, 



Aj{x) c|(£)c|(x -I- e^)- t<->i 



h.c. 



(A16) 



where the sum includes al l latt ice sites and j — 1,2. To 
ease comparison with Sec. [I A, we have returned now to 
a Hamiltonian formalism, replacing Cq, by c]j . We are now 
in a position to take the continuum limit. In this case, it 
suffices to neglect the massive 9 mode, and simply take 
Ai = — A2 — Aoe*"^. For agreement with Sec. |II A , we 



define = 2%/2Ao, or Ai = -A2 = A/2V2. In addi- 
tion, we take the continuum limit of the electron fields, 
using the decompositions 



*1 



112' 



-*12lH 



\x-y _ 



•.x+y _ 



-*^12 



-\T/ j^~y 

^222* ' 

-\T/ j^-y 
^221' ' 



^2 = J I cos(0i(x)-(/.2(a-)) + cos(02(f)-</'i(f-ei)) K^;^™*, with 



and the herniitian conjugates of these equations. In- 
serting these into Eq. Al£, gradient-expanding the 
fields, and rotating 45 degrees to a:; — y coordinates along 
the (7r,7r) and (— 7r,7r) directions, one obtains Hint — 



-|-cos((/)i(x-ei)-02(^-e2)) + cos((?!)2(^-e2)-(?I>i(:?)) |, 

(A13) 

with 7 = (Ao)^7. We assume 7 > 0, in which case this in- 
teraction favors a relative phase difference of 4'i~4'2 = ""j 
i.e. d^2_y2 order. 

We now turn to the fluctuations around the uniform 
dSC state. To do so, we let 0i = ip, (j)2 = + 9 + n. 
The phase ip is prec isely the order parameter phase in- 
troduced in Sec. HA. The other variable 9 represents an- 



1 <-> 2, a; <^ y). 



(A17) 



other branch of massive fluctuations around the d-wave 
state. We can thus assume 9^1, and that p is slowly 
varying, i.e. dj(p ^ 1. Under this assumption, we can 
take t he con tinuum limit and replace the positional sum 
in Eq. A13 by an integration. This gives H2 — J cPx'H2, 
with 



^2 = 7<! TT + li9jV>f - \dxpdyp + - dy)Lp 



This form is identical to the A term in Eq. 2.10| when 
the order parameter A is constant, but the symmetric 
placement of derivatives is important in the presence of 
phase gradients. In particular, now let A = v^e^'f and 
integrate by parts to transfer the derivative in the second 
term from the ^'^ to the A^ combination. Then, using 
the operator identity 

i {e'^idy + idye'^P) = e'^P'^idye'^P'^, (A18) 
one obtains 

n,nt = *I[^^Ar^e'^^/2(i9,)e-^/2]vE-i + (1 ^ 2;a; ^ y). 



s=± 



(A19) 



Eq. A19 is identical to the symmetrized form of the 



(A14) phase-quasiparticle interaction hypothesized in Eq. 2. If; 
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APPENDIX B: RENORMALIZATION GROUP 
ANALYSIS 

In this appendix, we present some details of the RG 
calculations for the transition at half-filling with parti- 
cle/hole symmetry between the nodal liquid phase and 
the antiferromagnet. As discussed in Section IV, this 
transition is described by the same field theory as the 
phase transition between the d-wave superconducting 
phase and the phase with co-existing antiferromagent and 
gapped d-wave order. A striking feature of this transi- 
tion, which we access perturbatively in an e = 4 — (d+ 1) 
expansion, is that it is Lorentz- invariant. Although the 
model possesses three independent velocities, the differ- 
ences between these scale to zero at the critical point. 

The full Lagrangian is given by. 



2 „__ivt2 -.(']Sf2\2 



(Bl) 



where we have suppressed the spin subscripts a, (3 on the 
<T Pauli matrices, and the particle/hole subscripts a, 6 on 
the T Pauli matrices. Here Kq — K and Kj = —v^K, for 
j = 1,2. Notice that this model has three independent 
velocities, vp,v^ and Vs- 

Now, we can rescale the N field to set if = 1, and 
rescale x to set Vs = 1- The Lagrangian can then be 
rewritten as: 



lidtNf 



- (;N • [ip'^T^aa^'ip^ + C.C.] 



(B2) 



where 



ai = 1, 

Vs 

02 = 1. 



Power counting about the Gaussian theory reveals that 
both u and are relevant in £> = 2+1 = 3 space-time di- 
mensions, but become marginal in Z) = 4. To implement 
a perturbative RG calculation thus requires continung 
the model above D = 3. This is a little tricky, due to the 
Pauli-matrix algebra. One approach is to dimensionally 
continue the loop integrals, but leave the 3-D Pauli-matrix 
algebra unchanged. This turns out to be equivalent to in- 
troducing into the Dirac equations an extra Pauli matrix, 
T^, multiplied by another (3rd) spatial dimension. Alter- 
natively, one could replace the two-dimensional Pauli ma- 
trices by the 4— dimensional 7-matrices, appropriate for 
4-D spinors. This latter procedure would be more correct 



if we were truly interested in the vicinity of four dimen- 
sions, but our choice probably makes more sense given 
that we are eventually concerned with e ^ 1. In any 
case, the difference is fairly trivial: factors of 2 would be 
replaced by factors of 4 in traces over the 4D 7-matrices. 
With our convention, we obtain the following one-loop 
flow equations: 



d''^ AA 2 

— = eu - 441*"^ 

dl 
dl 

dvN _ 1 

— - = - TAT = 2 

dl V 



32Au + 32A2 



20u- 16A)rjv 



(B4) 



where X = g^. N and ip have anomalous dimensions 
8A and 3A, respectively. These flow equtions have a fixed 
point at 0(e): A = e/20, u = e^/37T/440, tn = 0, = 0. 
This fixed point has the following interesting features: 

(1) 'Relativistic Invariance.' Since the velocity differ- 
ences scale to zero according to ( B4 ) , all physically mea- 
surable quantities are a function of x^+y'^—t'^ in the units 
which we have chosen, or, upon restoring the velocities, 
{x^ + y^)-v^,e. 

(2) An antiferromagnetic correlation length which di- 
verges as 



(B5) 



as the transition is approached. 

(3) Critical correlation functions with the following 
power-law decays 



{N,{x,t)N, {0,0)) 
(V't(x,t)V'(0,0)) 



|f2 _i2|i-16A* 
1 

|f2 _^2|1-6A- 



(B6) 



As a check on the reliability of tlie e-expansion, we 
consider the following related modeled: 



(B3) c= i^dt 



T^idx 
1, 



+ l^mt - ^la,*!' - ^^1*1' + <?[^VV + C.C.], (B7) 

where -0 is now a single two-component spinor, and $ is a 
complex field. Using the same e-expansion procedure, we 
find a fixed point at X* ^ u* = e/6; at this fixed point, 
the fields $ and ij) have anomalous dimension e/6. It is 
a remarkable and fortunate fact that this model exhibits 
= 1 supersymmetry. As a result, the existence of a 
fixed point is guaranteed and the scaling dimensions of tp 
and $ can be determined exactly, in agreement with the 
e expansion. 
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